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1 Introduction 



Great strides have been made recently in automatic summation of series, particularly 
liypergeometric series. A source for this is jPWZ96j . which includes a historical 
development of the problem as well as a fine exposition of the recent and seminal work 
of the three authors. A cornerstone of the automation of hypergeometric summation 
is Gosper's algorithm. In )MG02aj . it is pointed out that "despite the fact that 
Gosper's algorithm is one of the most important achievements in computer algebra, 
to date there are no results concerning the average running time of the algorithm." 
In that same work, Milenkovic and Compton undertake an analysis of the run time 
under various stipulated probabilistic models for the inputs. 

To describe the results of MCI, let / and g be polynomials, let := f{k)/g{k), 
and let t„ := 11^=1 '^fc- The series {tn} and its partial sums are known as hypergeomet- 
ric, and Tk is called the hypergeometric ratio. The purpose of Gosper's algorithm is 
to find a closed form expression for the partial sum Sn '■= '^n=i ^n- Its input is often 
specified as the rational function rk in factored form. Roots of / and g differing by 
integers play a crucial role in the algorithm. Milenkovic and Compton observe that 
not much generality is lost in assuming the roots have been classified according to 
their remainders modulo 1, and that the problem has been restricted to one of these 
moduli classes. In other words, they assume that / and g have integer roots. They go 
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on to stipulate joint probability distributions for / and which have as parameters 
an a priori bound on the location of the roots. Specifically, they assume that 

m 

i=i 

m 

9{k) = \{{k-3r^ (1.1) 
so all roots of / and g lie in [m] :={!,..., m}. 

The "uniform R model" considered by Milenkovic and Compton may be described 
as follows. For each of / and g, a sequence of n IID uniform picks is made from [m] . 
This gives the multisets of roots for / and g. In other words, for j e [m], the 
random variables Aj and Bj count how many of n IID picks from [m] are equal to j. 
Milenkovic and Compton point out that this is an urn model of Maxwell- Boltzman 
type. A key to their analysis is the representation of the variables {Aj} as a set 
of m IID picks from a Poisson distribution (of any mean), conditioned to sum to 
n. A second distribution of roots they analyze is the "multi-set R model" , in which 
the multisets of roots (equivalently the sequences (^j)^i and (-Bj)^i) are chosen 
uniformly from all multisets (equivalently all sequences of m nonnegative integers 
summing to n). This is a Bose- Einstein urn model, and is equivalent to conditioning 
two IID sequences of geometric random variables (with any mean) both to sum to n. 
They also discuss two models, the "Uniform T model" and the "Multiset T model" , 
in which the roots of the numerator and denominator of the partial product terms 
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are directly modeled by the two respective urn models; these models are not addressed 
in this paper. 

Milenkovic and Compton give a partial average case analysis, meaning that they 
focus on a few quantities which are highly determinative of the run time and give 
average case analyses of these. The most important such quantity is the degree of the 
cancellation polynomial, c{k). This is defined as the minimal polynomial c for which 
we may write 

m ^ aik) cik + l) 

g{k) b{k) c{k) ^ ■ ^ 

and also satisfying 

GCD(a(a;), h{x — h)) = 1 for all nonnegative integers/i . (1-3) 

The determination of this polynomial is Step 2 in the version of Gosper's algorithm 
described in |Wis03j . which is distilled from |PWZ96j . Milenkovic and Compton 



obtain the following results (they use in place of the m in this paper). The draft of 
their manuscript cited here is a very preliminary version which the authors have kindly 
provided. Consequently, only results independently proved in the present paper are 
quoted here, though in fact the manuscript |MC02aj obtains explicit constants for the 
asymptotic expectations. 

Theorem 1 (Milenkovic and Compton (2002) Theorems 20 and 21) In the 
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uniform R model, ifn/m X, the expected degree of c{k) is asymptotic to 

In the multiset R model, they find that when A is sufficiently large, the expected degree 
of c is asymptotically 

The method of |M(]n2aj is to compute transforms (generating functions) for the 
unconditioned distributions, in which the variables Aj and Bj are independent Pois- 
sons or geometries, and then de-Poissonize, according to machinery they developed 
M(]n2bj . 
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Analytic de-poissonization may be technically somewhat involved; see for exam- 
ple |JS99j . The present paper also relies on the representation of the stipulated 
distributions as IID conditioned on a fixed sum. After that, however, the method 
herein is purely probabilistic, relying on limit theory for the random walk whose in- 
crements are Aj — Bj. Theorem El below, whose proof is a straightforward application 
of random walk limit theory, shows that as n, m oo with n/m = A + o(A^^/^), the 
expected degree of c is 

(c + o(l))m3/2 ^^4) 
where c is a certain expectation taken with respect to the Brownian bridge. For the 



uniform and multiset R models, the constant c is calculated respectively as 



Cunif = ; (1.5) 

lb 

^V2A(A + 1) 

Cmulti = 77 ■ (1-6) 

io 

The authors of M(^n2al are aware of the random walk representation, but it appears 
that they use this only via analytic transforms, and not via any scaling limits of the 
random walk paths. 



2 Definitions and results 



Let F be a distribution on the nonnegative integers with mean A and variance V < oo. 
We assume throughout that the GCD of the support of F is 1, as is true for Poisson 
distributions and geometric distributions of any mean. Let P_F,m be the probability 
measure on S := Z^"* making the coordinates IID with common distribution F. De- 
note the first m coordinates hj Ai, . . . , Am and the last m coordinates hj Bi, . . . , B^- 
Let QF,m;n bc the result of conditioning P_F,m so that 

m m 

Associated to each uj E S are the polynomials / and g defined by and the 

associated hypergeometric series with ratio = f{k)/g{k). Define a function /3 : 
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S Z+ by letting f3{uj) be the minimal degree of a polynomial c satisfying (|1.2p for 
the polynomials / and g. The next result, proved at the end of the section, provides 
an alternative expression for (3, which is essentially the random walk representation 
in |MC;()2al Theorem 1]. 

Define Xj := Aj — Bj and Sk '■= Y2'j=i-^j^ with the convention that 5*0 := 0. 
Define Mk := minjS'j : < j < A;} and define Mk := min{Sj : k < j < m}. Let 
r := min{/c : Sj > Sk Vj > k} denote the time of the first minimum of the process 
{Sj : < j < m}. For each j, define Yj := Sj — Mj and Yj := Sj — Mj. Define 
Ij := Yjlj^r + Yjljyr- At the end of this section we will prove: 

Lemma 2 

m—1 T m—1 

j = l j = l j=T+l 

Under the measures PF,m and QF,m;n, the quantities Aj, Bj, Xj, Yj, (3 and so forth 
become random variables. In order to state the main results of this paper, some 
definitions are required that mirror the definitions of these quantities but on the 
space of continuous limits. 

Let VL be the space of CADLAG paths on [0, 1], with filtration {J^t}- For G 
define the minimum process M* by M*(uj){t) = inf{co'(s) : s < t}. Define M* to be 
the right to left minimum process M*{uj){t) = inf{ci;(s) : s > t}. Define the process 

6 



Y toheu-M* and Y toheu- M*. Define r = mi{t : u;(s) > u;(t) Vs > t} to be 
the time of tlie first minimum and let / be tlie process Ylt<T + Ylt>T- Finally, we 
define P*{u;) to be the quantity I{u){t) dt. 

Let P2V be the law of a centered continuous Gaussian process with covariances 

Eu;{s)u;{t) = 2V{s At) . 

In other words, P2V is a Brownian motion with amplitude \/2V. Let Q2y be the law 
of a centered continuous Gaussian process with covariances 

Euj{s)uj{t) = 2V{{sAt)-st). (2.7) 

In other words, Q2V is a Brownian bridge of amplitude y/2V. Recall that the mean 
and variance of F are denoted by A and V. The first main result of the paper, proved 
in the next section, is: 

Theorem 3 Let m — > 00 with n — Am + o(m^/^). Then 

J P d^F^m = {Ki + o(l))m^/2 (2.8) 

and 

J P dQF,m;n = (i^2 + o{l))m^/^ (2.9) 

with 




(2.10) 
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K2:= j P*dQ2V. (2.11) 
The integrals ()2.10|) and ()2.11|) may be evaluated, leading to quantitative versions: 



Theorem 4 Consider the measure PF,m, under which {Aj,Bj : 1 < j < m} are IID 
and have common variance V . Then as m oo, 

E(3 = {1 + o{l))l\l^m^/' . 

3 V TT 

Theorem 5 Suppose that N{m) is an integer satisfying n{m)/m = A + o(m^^/^) as 
m — >• oo. Consider the measure ClF,m;n, under which {Aj,Bj : 1 < j < m} have the 
distribution of IID picks from F conditioned on Y^jLi ~ Xljli ~ ^ ■ Then as 
m oo, 

E/3=(l + o(l))^v^m3/2. 
Id 

Remarks: 

(1) Theorem El emerges without much difficulty from the convergence of 
the random walks to Brownian paths. Thus not only does m~^l'^ times 
the PF,m expectation of (5 converge to E2y/5*, but the P_F,m distribution 
of mr'^l'^(5 converges to the distribution of 
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(2) Although it is relatively easy and is superseded by the two quantitative 
results, TheoremElis worth stating separately for the following reason: the 
computations in Theorems |3] and El are a little tricky, and it is instructive 
to see that the form of the result does not depend on calculations which 
are not transparent. 

(3) The asymptotics obtained by |MC02aj for both the uniform and mul- 
tiset models pertain to Theorem Theorem |3] corresponds to two models 
discussed in |MC02aj but not quoted here, where the roots of / and g 
are assumed to be unconditioned picks from the Poisson (respectively ge- 
ometric) distributions, hence not necessarily equinumerous. 

The rest of the organization of this paper is as follows. The functional central 
limit arguments are spelled out in Section El The values of the constants Ki and K2 
are then computed in Section 01 In the remainder of this section we prove Lemma |21 
We begin with an intermediate representation. 

For each j between 1 and m, place Aj red balls and Bj blue balls in an urn marked 
j; this will be called "position j" or "time j". An admissible matching of the balls is 
a set of pairs of balls such that 

(i) each pair contains one red ball and one blue ball; say that these two balls are 
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matched, and any ball not in the union of pairs is called unmatched; 

(ii) the pairs are pairwise disjoint; 

(iii) if a red ball in position i is matched with a blue ball in position j, then i < j; 

(iv) if there is an unmatched red ball in position i and an unmatched blue ball in 
position j then i > j (that is, among unmatched balls, all red balls sit to the 
right of all blue balls). 

For each admissible matching C,, define the weight w{C,) to be the sum over all pairs 
in the matching of j — i, where j is the position of the blue ball and i is the position 
of the red ball. 



Proposition 6 

f5 = min{w{C,) : C, is an admissible matching} . 



Proof: For every matched pair in positions i and j, we have 

X — i X ~ s e{k + 1) 

X — j ^J: X — {s + 1) e{k) 

where e{x) = YliZli^ ~ s)- Thus any admissible matching ^ of weight w yields 
a solution H{^) to p.2|) where c has degree w. Conversely, let c solve ()1.2|) and 
have degree w. If w = then the empty matching is admissible. Assume now, for 
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induction, that w > and that for any /', g' represented by red and blue balls in 
urns, and any solution c' to ()1.2|) of degree less than w, there is an admissible ^ for 
that urn problem with c' = H{^). For some j which is a root of c, let c' denote c 
with the linear factor {x — j) removed. By induction, there is a matching ^' of weight 
w — 1, with H{^') = c', admissible for the urn problem gotten from the original one by 
adding a blue ball at position j and a red ball at j + 1. Both of these new balls must 
be matched in say in pairs of positions (s, j) and (j + l,t), whence replacing these 
with the pair {s,t) produces a ^ admissible for the original problem with H{^) = c. 
□ 

Proof of Lemma |21 For any admissible matching ^, define dj{^) to be the number 
of red balls in positions 1, . . . , j matched with blue balls in positions j ' + 1, . . . ,m. 
Elementarily, summing by parts, 

w{o= J2 Ei = E^*(^)- (2.12) 

pairs of J<i<« t 

positions 

Fix any j < r and let = argminj^j : < i < j} so that Sj^ = Mj. Then 
in positions + 1, . . . , j there are a total of Yj more red balls than blue balls. In 
positions j + 1, . . . , r, there is an excess of 5*^ — Sr > Sj — Mj = Yj blue balls over 
red balls. In an admissible matching, either every red ball at a position at most j 
is matched or every blue ball at a position at least j ' + 1 is matched. It follows that 
either at least Yj red balls from positions + 1, ■ ■ - j are matched to blue balls in 
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positions beyond j, or at least Yj blue balls from positions j + 1, . . . ,r are matched 
with red ball in positions up to j. In either case, dj{^) > Yj. 

Similarly, fix any j > t and let j* > j satisfying Sj* — Mj. There is an excess of 
Yj-i blue balls in positions j, . . . ,j* — 1. There is an excess of red balls in positions 
T + 1, . . . , J — 1 of Sj^i — St- > Yj_i. Reasoning as before, one sees that > l^-i. 

Summing over j now gives 

m— 1 T m—l 

J2djim)>J2^j + J2^j. (2.13) 

j = l j = l j=T 

Note that Y^ — Yr — 0, so this is equal to Yl]^=i ^j- Minimizing over ^ then gives half 
of the conclusion of the theorem: /? > Yl^=i 

To prove the other half, we produce an admissible matching ^ with dj{^) — Ij for 
sill < j < m — l. In particular, dr = 0, so no ball in a position at most j is matched 
with a ball in a position beyond j, and ^ may be decomposed into a matching on 
balls in urns 1, . . . ,j and another on balls in urns j + 1, . . . ,m. We construct these 
separately. An algorithm for the first is as follows. 

Initiahze j :— 1. Initialize a LIFO stack. Pull red balls out of urn j and 
place them on the stack until there are no more red balls in urn j. Pull 
blue balls out of urn j: while the stack is non-empty, match each blue ball 
with the top element of the stack; once the stack is empty, label each new 
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blue ball "unmatched" and discard it. When urn j is empty, increment j 
and execute the loop until finished with the step j = t. 

It is easy to see that all red balls in positions up to r will be matched, since every 
time a red ball goes on the stack, there are more blues and reds to follow by time r 
and no blue will be discarded until that red ball is matched. Inductively, it is easy to 
check that: 

• the stack size after step j is Yj; 

• the change in stack size from time j — 1 to j is max{Xj, —Yj}; 

• the total number of balls discarded through time j is —Mj. 

From this one sees that dj{^) is equal to the stack size after time j, and is therefore 
equal toYj. A stack algorithm dual to this works in the case j > r, working backward 
from time m to r, stacking blue balls and matching or discarding reds. It constructs 
the other half of C, so that dj{C,) > Yj for all j > r. This completes the proof of 
Lemma |21 □ 
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3 Proof of Theorem |31 



Let T be the topology on Q generated by the sup norm := supo<t<i |^^(^)|- The 
following lemma is necessary only because r is not a continuous function. 

Lemma 7 (3* is a continuous function on Vt with respect to T . 

Proof: Suppose sup^ — uj2{s))\ < e. Fix any t. It is clear that \M{t,uJi) — 

M{t,uj2)\ < e and likewise for M, hence these are continuous. If r(ci;i), r(co'2) > t, it 
follows immediately also that \I(t,uJi) — I(t,uj2)\ < 2e. Likewise, if t{ui),t{u!2) > t, 
it follows that |/(t,cui) - /(t.cus)! < 2e. 

Suppose now that Ti := r(co'i) < t < T2 := t{uj2)- Then 

M{t, ui) > uiin) > uj2{n) - e > M{t, U2) - e 

since ti is one of the times over which the inf defining M{t, 002) is taken. Similarly, 

M(t, U2) > uj2{t2) > cji(r2) - e > M{t, Ui) - e . 

It follows that 

M{t,tui) - M{t,uj2) <e. (3.14) 

Together with \u!i{t) — uJ2{t)\ < e, this shows that \I(t,u!i) — I(t,U2)\ < 2e. A similar 
argument shows this in the case that T2 < t < ti. This establishes continuity of /, 
with continuity of D following by integration. □ 
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Proof of (j2 . 8j) : Recall the definition of the partial sums Sj on S and define a map 



are 


evident: 








T* O K 


= m~ 


-V; 


M* 




= m~ 




M* 




= m~ 




Y* 




= m~ 


-^'^Y{Lo{[mt\)) 


Y* 


[K{uj){t)) 


= m" 




r 




= m' 






[3* OK 


= m~ 





Let P(m) denote the image under k of P^.m- The functional central limit theorem 
says that the laws under P(m) of n converge weakly as m — oo to the measure P2y- 
See |Bil86t Theorem 37.8] for a proof when p > 4 or |Dur961 Theorem 6.3 of Chapter 7] 



for a general proof using Skorohod embedding. This and Lemma [7| would complete 
the proof of (j2.8p if (3* were bounded. Since (3* is not bounded, we may define for 
each L > 0, 

/i:=sgn(r)(|r| AL) . 
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We may then conclude that the expectation with respect to P(m) of the bounded 
continuous function i3{l) '■= J I{l) converges as m ^ oo to its expectation with 
respect to P2i/. 

Lemma 8 Lei ^max := supo<t<i Then 

hm sup E(m)^maa;l/^^g2>L = 0- 

Assuming this for the moment, we observe that P(m) is symmetric so the same 
holds with inf in place of sup, and thus 

lim sup E(^)Bspanlg3 ^^^^ = (3.15) 

where fispan = supjCi;(t) — 'm.ituj{t). Since 

— (3(^L){ijj)\ < Bspan := supu;(t) — infti;(t) 

t * 

and 13* = on the event {Bspan < L}, the inequality (j2.8p follows from ()3.15p and 
the convergence of Fj(^rn)Pp{L) to E2vP{l)- 

Proof of Lemma |H1 By the maximum inequality ( |Dur96| Theorem 4.4.3]), 



E(^)Bmax2 < AEu{lf 



-'m I 



4Ei.,„(m-i/2^„ 
AV, (3.16) 
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and hence 



E(m)Bmaxlj3j^g^jj^^ < L ^E(m)Bmax^ < ^ (3.17) 



for all m, proving the lemma. □ 



Proof of (12. 9^ : Recall that Q,2V denotes the law on of a Brownian bridge of 
amplitude \/2V . The proof proceeds analogously to the proof of ()2.8|1 . In place of the 
standard functional central limit theorem is a well known result that may be found, 
among other places, in |Pit02| equation (6) of Section 0.4] (refer to |DK63j for the 
proof). 

Lemma 9 (Conditional Functional CLT) Letn{m) ^ oo asm —> oo with n{m) / m 
\ + o{m~^^'^). Recall that Q,F,m;n is the measure on (Z"*")^™ whose coordinates have the 
distribution of IID draws from F conditioned on Xljli — XljLi — ^'^'^ recall 
the aperiodicity assumption on F. Let Q(m) denote the image under k of Q,F,m,;n- 
Then the Q(m) lO'W of {Sj : 1 < j < m} converges weakly to Q2V- D 

All that remains is to show the analogue of Lemma |H1 with Q{m) in place of P(m). 
As in (|3.17p . this will follow once we have shown a uniform bound on the Q(m) second 
moment of Bmax^ analogous to (|3.16|) . This in turn follows immediately from the 
inequality 

QF,m;n(fimaX > L) < CP F,mi^'^8.X > L) (3.18) 
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where hereafter C may change from equation to equation but will depend only on F. 

To prove ()3.18|) . let G be the event that Yl^=i = Xljli = ^- the local central 

limit theorem, Pp^rn{G) is asymptotic to C/n, hence, using time-reversal symmetry 

of the path {Sj} under PF,m, 

, PF,m(Bmax > L; G) 
Qj7,m;„(fimax > L) = — 

< CnPir,m(Bmax > L; G) 

= 2Cn^ ^ PF,m{iimax = j = St;G) . (3.19) 

j>L t<m/2 

It then suffices to show 

PF,„(Bmax = j = St; G) < Cn-^P^,„(Bmax = j = St) (3.20) 
since then resumming fl3.19j) proves ()3.18|1 . 

Let / = [3m/4j, let J-'i be the cx-field generated by {Ai,Bi : i < 1} and let H be 
the event that St = maxj^j : i < I}. The local central limit theorem |Dur96t (5.2) in 
Chapter 2] gives PF,m{G \ J-'i) < Cn^^, whence 

PF,m{^max = j = St;G) < PF^St = j; H nG) 

< PF,m{St=j;H)PF,m{G\J'l) 

< Cn-'PF,m{St=j;H) 

Conditioning again on JF;, we see that PF,m{St = j; H) < Cn^^PF,m{St = j = Bmax), 
which establishes ()3.2()j) . hence ()3.18j) and the theorem. □ 



4 Evaluation of the constants 



Let 

Di{iu) := [ {uj{t) - M{t))dt. 
Jo 

By symmetry, / Di rfPay = (1/2) / /J* dP2v and similarly, / Di dQay = (1/2) / (3* dQay- 
Thus it suffices to compute expectations of Di. 

Proof of TheoremEI For a process with law P2V, the process Y{t) = uj{t)-M{t) 
is well known to have the same distribution as the law under P2V of the reflected 
Brownian motion {|co'(t)| : < t < 1} (see, e.g., |Kal02t Prop. 13.13]). Clearly, the 



map uj ^ uj ~ M has the property that t is a left-to-right minimum for uj if and only 
if t is a zero of u; — M. The last left-to-right minimum is the global minimum, which 
occurs at the last zero of a; — M. Because the process u; — M is distributed as 
the location of the last zero oi u — M is distributed as the last zero of uo. This has 
an arc-sine density 7r^^(x(l — x))^^/^ dx |Dur96t Example 4.4]. The following lemma 
writing Brownian motion as a mixture of bridges up to the last zero follows directly 
from the strong Markov property, scaling, and the fact that the bridge is a Brownian 
motion conditioned to return to zero: 

Lemma 10 Let L = L[lo) = sup{t < 1 : u{t) = 0} be the last zero of Brownian 
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motion and let g : Q ^ depend only on co'|[o,l((j)] ■ Then 

Jn Jo 7rA/t(l — t) Jn 

where the inner integral is on J^t and Q,2V ^-^ ^^''^ -^t ^/ ^ bridge of amplitude 
^/2V on [0,t], that is, a centered Gaussian process with covariance 



Euj[u)uj[v) = u Av . 



□ 



Using this, c may be evaluated as follows. By definition, by |Kal021 Prop. 13.13], 
and lastly by Lemma ITUl applied to the integral up to L of |ti;|, we have 

[ DidP2viuj) = [ [ iujis)-M{Lu){s))dsdP2V 
Jn Jn Jo 

= \uj{s)\dsdP2v 
Jn Jo 



^ dt f 



7T\/t{l - t) Jn "^^^JO 

From the covariance structure of Q2V', we see this law makes uj{s) is a centered 
Gaussian with variance s(l — s/t). The expected absolute value of a A^(0, a) random 
variable is ^/laji:. Switching the order of the two inner integrals, we may then write 

/ ArfP^vM = t ^^== f ds I \uj{s)\dd%{uj) 
Jn Jo T^\/t[l - tj Jo Jn 



dQiy / \uj{s) \ ds 



dt r ds lAV 



n^t{l - 1) Jo -Ky^til-t) V TT \ \ t 
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The evaluation is now straightforward integration. Substitute s = tu and ds = tdu 

to get 




Doubhng yields / /?* dP2v and finishes the proof of Theorem HI that is, Ki = 
(2/3)^^7^^ 0.376v^. □ 

Proof of Theorem |S1 If the distribution of uj{t) — m{t) for a Brownian bridge 
were explicitly known in a usable form, the computation of the Q-expectation would 
be analogous to the P-expectation of D. In the absence of such a representation, the 
second computation ignores the representation of the law of — M as that of \uj\ and 
proceeds as follows. 

The counterpart to LemmalTUlis the following joint density for the pair {uj{t), M{uj){t)) 
under Qi conditioned on r > t. 

Lemma 11 For fixed t G (0,1), define the positive function fn on the set R : = 
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{{x,y) : y >0x > —y} as follows: 

Then ffj is a conditional density for {uj{t), —M{uj){t)) under Qi conditioned onr > t. 



Proof: Begin with the computation of a density for (u;, M) under the standard 
Brownian measure P. By the reflection principle, using P" to denote standard Brow- 
nian motion starting at a, one has 

P^{uj{t) e [x,x + dx],M{t) < -y) = P''^y{uj{t) G [x,x + dx]) 



2TTt 



Differentiating with respect to y yields 



pO(^(i) e [x,x + dx],Mit) e [-y,-y + dy]) = ii^±M^_Le-{-+2,)V(2t)^^^^ 

(4.21) 

On R. 

Next, compute the joint density 

P{uj{t) G [x,x + dx],M{t) G [y,y + dy],uj{l) G [0,dz],T > t) . 
To do this, according to the Markov property, one must multiply ()4.21|) by 



P^(a;(l-t) G [Q,dz] min uj{s) < -y) 

\)<s<l—t 
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By the reflection principle, this last factor is equal to 



p-^y-(uj{l -t)e [0, dz]) = ^-l-e-(^-+2^)'/(2(i-t)) ^ 

and multiplying and simplifying 1/t + 1/(1 — t) to l/(t(l — t)) in the exponent gives 
a joint density of 

x + 2y ^^_(,+2,)2/(2,(i_,)) _ ^^ 22) 



TT 



A change of variables simplifies the computation a little. Let u = x + 2y and 
f = (2x — ?/)/5, so that c/f = c??/ and x = {u + lOf )/5, y = {2u — 5f )/5. The 
region R is transformed into the region R' := {m > 0, (— 3/5)m < f < (2/5)ti. Now 
rewrite the density ()4.22p as 



« e--y(ni'^))dzdudv. (4.23) 



The conditional density of {u,v) given uj{l) G [0,6/2;] and r < t, is given by 
normalizing this. One must divide by the integral of (j4.23|) over R', computed by a 
simple linear change of variables u = r-^t(l — t) in the third line: 

r-oo /.(2/5)« 



J-{3/5)u TT>yt^{l - t) 



U 



2 



e 



TX^t\l-t) 

t(l-t) e-'''^Jt{l-t)drdz 

— 1= dz 
V2^ 
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using the fact that r^e ''^''^ dr = ^/ttJ2. Dividing, 



and plugging in u = x + 2y proves the lemma. □ 

Proof of Theorem O continued: Let G denote the CDF for the time at which 
a Brownian bridge on [0, 1] reaches its minimum. Set the amplitude 2V = 1 for 
convenience, and note that co'(t) — M(t) = x — {—y) = {3u+5v) /5. Then bv LemmalTH 
we have 

[ DidQi = If {uj{t) ~ M{t))lr>tdtdQli{uj) 
Jn Jn Jo 

= [ dt{l-G{t)) /"(cj(t)-M(cu)(t)) rf(Qi|r>t)(cj) 
Jo Jn 

dt{l-G{t)) I — ^ — fH{x{u,v),y{u,v))dudv. (4.24) 



JR' 



The integral over R' may be computed by substituting u = r^t{l — t) as before 
to get 



iJ' 5 ' y 7ri:3(l - tf Jq 7(-3/5)n 5 



oo 



3 o 1 



7rt3(l -t)3 ^5 10 ' 
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Finally, one may evaluate the integral in (j4.24j) without knowing the exact distri- 
bution function G. By symmetry, we know G{t) + 0(1 — t) = 1. The integral over R' 
is symmetric under t ^ 1 — t. Thus 

f D,dQ, = f\l-G{t))\^/t{l^)dt 
Jn Jo 2 

71 

32' 

Now doubling to get D and multiplying by the amplitude of V2V leads to 

Jn 16 
which establishes the value of K2 and finishes the proof of Theorem 
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